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With the use of the relativistic string Hamiltonian we study the universal static 
potential V s t(r) and the force, which are fully determined by two fundamental pa- 
rameters: the string tension a = 0.18 ±0.02 GeV 2 and the QCD constants A^jg(nj), 
taken from pQCD, while the infrared (IR) regulator Mb is expressed via the string 
tension. The vector couplings ay(r) in the static potential and ap(r) in the static 
force, as well as the characteristic scales, r\{nf = 3) and ro(nj = 3), are calculated 
and compared to lattice data. The result roAy^(nj = 3) = 0.77±0.03 is obtained for 
Mb = (1.15 ± 0.02) GeV, while in bottomonium better agreement with experiment 
is reached for A^jg(n/ = 3) = (325 ± 15) MeV and the frozen value ay = 0.57±0.02. 
The mass splittings M(LD) - M(1P) and M(2P) - M(1P) are shown to be sen- 
sitive to the IR regulator used. The masses M(1 3 D 3 ) = 10170(2) MeV and 
M(l 3 L>i) = 10154(3) MeV are predicted. 



I. INTRODUCTION 

The Hamiltonian formalism may be considered as a powerful tool to study such hadron 
properties as meson spectroscopy, including high excitations, hyperfine and fine structure 
splittings of different meson multiplets, leptonic widths, and radiative and strong meson de- 
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cays. For decades, different phenomenological Hamiltonians were used in constituent quark 
models, and some of them were rather successful in predictions of meson properties for low- 
lying states . However, in such models the quark- ant iquark potentials contain a large 
number of arbitrary parameters like constituent quark masses, variable values of the string 
tension and the QCD constant A, as well as a overall additive fitting constant. Meanwhile, 
the relativistic string Hamiltonian (RSH) H^, which was derived from the gauge-invariant 
meson Green's function with the use of the QCD Lagrangian |6], contains only fundamental 
parameters: the quark current masses, the string tension a fixed by the slope of the Regge 
trajectories for light mesons, and the QCD constants A(nf), taken from perturbative QCD 
(pQCD). It is important that in RSH the spin-independent static qq potential V st (r) is appli- 
cable for the quarks with arbitrary masses (including m q = 0). It is defined via the vacuum 
average over the Wilson loop (W(C)) and has to be a universal. In Refs. 6|- 9] the only 
approximation made is that (W(C)) is taken in the form of the minimal area law, which 
appears to be a good approximation for the separations r > T g ~ 0.15 fm, where T g is the 
vacuum correlation length. Later the behavior of the strong coupling in the IR region was 
studied in Ref. jioj], where it was shown that the IR regulator is not a new parameter but 
expressed via the string tension. 

It was shown that the nonperturbative (NP) part of the static potential has a linear 
behavior beginning at r > 0.2 fm, while at short distances, r < 0.15 fm, the NP potential 
is proportional to r 2 Such a deviation from linear behavior, in a very narrow region, 
gives a small effect for all mesons, with the exception of T(1S), which has a small size, 
R ~ 0.20 fm, and for T(1S) some corrections to the confining potential should be taken into 
account. 

This universal quark-antiquark potential is additive, containing the linear confining part 
and the gluon-exchange (GE) part: 

V st (r) = ar + V GE (r), (1) 



which is confirmed by the Casimir scaling, studied analytically [12(] and numerically on the 
lattice 131 ] . Here the string tension is fixed by the slope of the Regge trajectories of the light 
mesons, being known to good accuracy, a = 0.180 ± 0.002 GeV 2 . 
In the gluon-exchange (GE) potential ([1]) 

w»w = CO 



3 



the vector coupling in coordinate space ay(r) , is defined through the vector coupling ay(q 2 ) 
in the momentum space, 

oo 

. 2 f sin(gr) 

a v{ r ) — - / dq o>v[q). (3) 

J q 
o 

For large q 2 there exists an important relation between ay(q 2 ) in the momentum space and 



the conventional a s (q 2 ) in the MS renormalization scheme [14J. In pQCD the cross sections 
and other observables are predicted in terms of this coupling. 

The coupling a s (q 2 ) is measured at different (large) energy scales q 2 and the values 
obtained are usually presented at a common energy scale, equal to the Z-boson mass, 
Mz = 91.188 GeV. From numerous experimental studies, like the hadronic widths of the Z° 
boson, the r-lepton decays, radiative T(liS') decays, jet production in the e + e~ annihilation, 



and the structure functions in deep inelastic scattering, the world average value of t 



15 



re strong 
As 



coupling is now determined with a good accuracy, a B (mz) = 0.1184 ± 0.0007 
a consequence, the QCD constant Aj^g(ny = 5) is also know with a good accuracy. Then, 
using the matching procedure at the quark mass thresholds, the other A^inj) for rif = 3, 4 
are calculated and the three- loop calculations give the following Aj^g(n/) [1 61 ] : 

A m (n f = 3) = (339 ± 10) MeV, 
A m (n f = 4) = (296 ± 10) MeV, 

A m (n f = 5) = (213 ± 8) MeV. (4) 

This result is also important for the vector coupling since the "vector" constants Ay (nf) 
are defined via Ay^(rif) They appear to be significantly larger, e.g. Ay (rif = 3) = 

(500 ± 15) MeV corresponds to the value A m (n f = 3) = (339 ± 10) MeV from Eq. fll) (see 
Sect. ED. 



The analysis of the Vnv,(r) shows that perturbative effects determine ay(r) only at very 
small distances r < 0.06 fm 17J and this result was confirmed by the lattice measurements 



of the static potential [18] . It means that in the GE potential NP effects become important, 
beginning from very short distances, and one needs to determine the vector coupling in the 
infrared (IR) region. 

For heavy quarkonia the importance of NP effects was understood already in 1975, just 
after the discovery of the charmed quark, when the Cornell group introduced the linear + 
Coulomb potential with rather large vector coupling, ay = constant = 0.39 [l], in which 
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the asymptotic freedom (AF) behavior was neglected. However, future studies have shown 
that the AF behavior of the vector coupling is very important, in particular for the wave 
functions (w.f.) and its derivatives at the origin 19], 20]. When the AF behavior is taken 
into account, even larger asymptotic values of ay were used in phenomenological models, 
e.g. the value ~ 0.60 in Refs. [3], 20J and such a large spread in these values occurs because 
of a different choice of the fitting parameters. 

On the fundamental level not many theoretical attempts were undertaken to determine 
the strong coupling in the IR region. On the phenomenological level a regularization of the 
strong coupling was suggested long ago, with the prescription to introduce the IR regulator 

2 /o 2 +M 2 \ 1 L 

into the logarithm In jg, changing it into In ( — j [21] . This IR regulator was interpreted 
as an effective two-gluon mass M2 g = 2m g . However, in QCD the appearance of the gluon 
mass is forbidden by gauge invariance and the meaning and the value of the IR regulator 
remained unsolved for many years. 



Recently a new result was obtained within background perturbation theory (BPTh) [101 ] . 
where indeed the type of the logarithm, as in Ref. 2l], was derived with the IR regulator 
(denoted as Mb), expressed through the string tension, so that Mb is not an additional 
parameter. This NP regulator has the meaning of the mass of the two-gluon system, con- 
nected by the fundamental string (white object) and its value is determined by the equation: 
M| = 2na. It gives M B ~ (1.06 ±0.11) GeV for a ~ 0.180 GeV 2 , the accuracy of the calcu- 
lations being determined by the accuracy of the WKB method used (~ 10%). In 22] the IR 
regulators were studied in so-called "massive" perturbative QCD, developed within Analytic 
Perturbation Theory. Their values obtained in two-loop calculations, also lie in the range 
(0.9 - 1.2) GeV. 

However, for a given Ay^(nf = 3) the variation of the regulator Mb in the range 1.0 — 
1.15 GeV gives rise to large differences in the asymptotic (frozen) value of the vector coupling, 
called a cr it in Ref. [sj]: acy(q = 0) = ay(r — >■ oo) = a ait (nf = 3). For example, taking the 
perturbative central value, Ajj^(nf = 3) = 339 MeV from Eq. (jlj) and the corresponding 
Ay {jif = 3) = 1.4753 Ajyfg(?2/ = 3) = 500 MeV, one obtains a C rit(2 — loop) equal to the large 
value 0.82 for Mb = 0.95 GeV and a smaller value 0.635 for the larger Mb = 1.15 GeV. 
Meanwhile, the critical value of the vector coupling is of special importance for the meson 
spectrum and therefore one needs to determine the IR regulator, as well as A^fg(n/ = 3), 
with great accuracy. 



Notice that smaller values of A v (% = 3) were used in Refs. and in pQCD Ay^(nf = 

3) = (292 ± 29) MeV, smaller than given in Eq. (@J), was calculated in Ref. |2J]. On the 
lattice the small value Aj^g(n/ = 0) = (245 ± 20) MeV was obtained in quenched QCD 25]. 

Our goal here is to study how the properties of the vector coupling and the bottomonium 
spectrum depend on the IR regulator and Ay(w/ = 3) used. We will show that the meson 
spectrum is governed by the ratio rf = ^ ( or rf = ^r~)- Then for large Ay(nf = 3) ~ 
500 MeV, as in pQCD, the value of the IR regulator M-q = (1.15 ± 0.02) GeV appears to be 



preferable, in agreement with the prediction from 10]. 



II. RELATIVISTIC STRING HAMILTONIAN 

The RSH was derived from the gauge- invariant meson Green's function, performing 
several steps (see . For a meson q±q 2 with the masses rrii and m 2 the RSH contains 

several terms, 

Hr = H + H S d + H str + H SE , (5) 

where the part i^sD refers to the spin-dependent potential, like hyperfme or fine-structure 
interactions; the term H str comes from the rotation of the string itself and determines the so- 
called string corrections for the states with / ^ 0, while Hq^, c omes from the NP self-energy 
contribution to the masses of the quark and the antiquark 2y]. All these terms appear to be 
much smaller than the unperturbed part H for all mesons, and therefore can be considered 
as a perturbation. The part Hq is derived in the form, 

#o = — + — + — + — + - J — + V st (r). (6) 

Here the variables Ui are the kinetic energy operators, which have to be deftermined from 
the extremum condition, |^ = 0, giving 



Ui = s J m f +p 2 (i = i )2 ), (7) 
and therefore Eq. (J2]) can be rewritten as 



H = y/ml + p 2 + + p 2 + V st (r). (8) 

This form of Hr is applied to heavy-light mesons, while for bottomonium it is very simplified. 
Among all mesons the bottomonium spectrum is of a special importance, since it has the 
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largest number of levels below the open flavor threshold. Altogether there are nine bb 
multiplets with Z = 0,1,2,3 and seven of them were already observed, and one can use this 
extensive information to test the static quark-antiquark potential. An additional piece of 
information on the coupling a s (/i) at different scales u may be obtained from studies of the 



hyperfine and fine structure effects in bottomonium 27 ]. 

By derivation, the quark (antiquark) mass rrii in the RSH is equal to the current quark 
(antiquark) mass, fhi{fhi) in the MS renormalization scheme, and therefore it is not a fitting 
parameter. In the case of a heavy quark one needs to take into account corrections pertur- 
bative in q R . i.e., to use the pole mass of a heavy quark, which is taken here to two-loop 



accuracy 



15j: 



= <1 + |^S)+ J (9) 



where £2 may be taken from 15j. For the b quark the pole mass can symbolically be written 
as mb(pole) = nib(nib)(l + 0.09 + 0.05), where the second and third terms come from the 
a s and a 2 corrections. In our calculations mfc(pole) = (4.81 ± 0.03) GeV is used, which 
corresponds to the conventional current mass m^mi,) = (4.22 ± 0.03) GeV. 

For bottomonium the calculated string and self-energy terms are very small, < 1 MeV, 
and therefore H R = H + H SD , where from Eq. (jSJ) H , 



H = 2^mt + pi + V st (r), (10) 

has a kinetic term, which formally coincides with that in the spinless Salpeter equation 
(SSE), often used in relativistic models with constituent quark masses. Such a coincidence 
between the kinetic terms in the SSE and the RSH, which was derived from first principles, 
possibly explains why the use of the SSE was successful in Refs. 

aaa. 

An important feature of the RSH is that it does not contain an overall additive (fitting) 
constant, which is usually present in models with constituent quark masses and also in the 



lattice static potential [24j. Notice that the presence of such a constant in the meson mass 
violates the linear behavior of the Regge trajectories for light mesons. On the contrary, 
with the use of linear Regge trajectories can be easily derived with the correct slope and 
intercept {9J (later we use the string tension a = 0.180 ± 0.002 GeV 2 , extracted from the 
slope of the Regge trajectories for light mesons). In heavy quarkonia low- lying states do not 
lie on linear Regge trajectories, because of strong GE contributions. 
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III. THE VECTOR COUPLING IN MOMENTUM SPACE 

First we consider the vector coupling acy(q) in momentum space, taking it in two-loop 
approximation, where the coupling does not depend on the renormalization scheme. For 
ay(q 2 ) we shall use the notation ctB^ 2 ), bearing in mind that it contains the IR regulator 
M B , determined as in BPTh Q: 



A)*B V A) *B 
(q 2 +M 2 )) 

where M B , entering the logarithm t B = In ^ — A a B , is not a new parameter but determined 

"1 A v 

via the string tension 



10|: 



Ml = 2na. (12) 

Here a is the string tension in the fundamental representation, a = (0.180 ± 0.002) GeV 2 . 
The accuracy of the relation (fT2"j) is determined by the accuracy of the WKB approximation 
used in Ref. fic3 |. which is estimated to be < 10%. Therefore 

M B = (1.06 ± 0.11) GeV. (13) 

The analysis of the bottomonium spectrum shows that the larger values, Mb = (1.15 — 
1.20) GeV, are preferable, if a large Ay(nj = 3) = (500 ± 15) MeV, corresponding to 
the pQCD value A^(nf = 3) = (339 ± 10) from Eq. (jlj), is taken, while for the smaller 
Mb = (1.05 ± 0.05) GeV and the same Ay one obtains a too large 2P — IP splitting and 
also a large b-quark pole mass, = 4.90 GeV. 

The vector coupling a^q 2 ) has the AF property at large momentum q 2 , where the "vec- 
tor" constant Ay(uf) is expressed through the conventional Aj^g(nj). This connection be- 
tween both As follows from the relation between the strong couplings (valid at large q 2 ), 
established in Ref. jl^ . 

avtf) . a , tf)(l + ^ a , (?) ) K _^M_. (14) 

Here a± = y — Y n f- For our purpose it is enough to use in Eq. (Ti~4"l) the one-loop approxi- 
mation for both couplings: Ob = ^ m jt with /3 = H — f 77 -/ an d the analogous expression 
for a s (q). From this equation it follows that In j^- = ^ln and therefore 

A v (n/) = A M g(n / ) exp • ( 15 ) 
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This relation gives 

A y (n f = 3) = 1.4753 A m (n f = 3), 

A v (n / = 4) = 1.4238 A m (n f = 4, 

A v (n/ = 5) = 1.3656 A M ( W/ = 5). (16) 

If one takes the perturbative Ajjg(rif) from Eq. (Hj), then the following values for the "vector" 
constants in pQCD are obtained: 

A v (n f = 5) = (291 ± 11) MeV, 
A v (n f = 4) = (421 ± 15) MeV, 

A v ( n/ = 3) = (500 ± 15) MeV. (17) 

For our purpose it is more convenient to make the matching procedure for the coupling 
«b(<7 2 ) in momentum space, not for a s (q 2 ). It is interesting to underline that in this case the 
calculated values of Av(ny) for nf — 4,5 practically coincide with their values in Eq. (flTl) . 
although now the IR regulator is taken into account. 

We give below two sets of Ay(rif) for two different values of Mb, equal to 1.15 and 1.0, 
namely, if M B = 1.15 GeV 

A v (n f = 5) = 310 MeV, A Y (n f = 4) = 429.6 MeV, A v (n f = 3) = 497.4 MeV, (18) 

and if M B = 1.00 GeV 

A Y (n f = 5) = 315 MeV, A y (n f = 4) = 435 MeV, A v (n f = 3) = 499.7 MeV. (19) 

One can see that the fitted values of Ay(nf) weakly depend on the IR regulator Mb, taken 
in the range (1.0 — 1.15) GeV, changing within ±5 MeV. (Here the matching was performed 
at the quark mass thresholds: q 5i = 4.20 GeV and g 43 = 1.50 GeV.) The difference between 
these two sets becomes manifest in the frozen value: a C rit(<? = 0,n/ = 3) = 0.630 for 
M B = 1.15 GeV and a crit = 0.819 for M B = 1.0 GeV. 

In Fig. CD we show the coupling a-Q^q 2 ) with Mb = 1.15 GeV and Ay{rif) from Eq. (fl~8l) . 
For a comparison we shall also use the set with a smaller A(nf = 3) = 465 MeV, 

A v (n f = 5) = 292 MeV, A v {n f = 4) = 406 MeV, A(n f = 3) = 465 MeV, (20) 

and a smaller a crit (g 2 = 0, rif = 3) = 0.5712. 
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FIG. 1: Compound a B (q) for M B = 1.15 GeV and Ay(n/) from Eq. CED and for M B = 1.00 GeV 
and Ay(rif) from Eq. (|19p . 



We give in Fig. [I] two curves for the compound «b(? 2 ) with the parameters taken from 
Eqs. dUD and (US]). 

The freezing phenomenon of the vector coupling was widely used in phenomenology 2j-[5] , 
where typical values a cr it ~ 0.54 — 0.60 were used. However, in the lattice static potential 
much smaller values a cr ; t = a(lat) were measured: a(lat) ~ 0.22 in quenched QCD (n/ = 0) 



281 ] and a(lat) ~ 0.30 in full QCD 29j. The reason for that discrepancy possibly comes 



so 



24 



rom an additional nor- 



30]. 



from lattice artifacts, present in the lattice GE potential, and a . 
malization condition, usually put on the lattice static potential 

Also in some lattice potentials saturation of the vector coupling takes place at very small 
distances, r ~ 0.2 fm 18| . 281 ] . while in our approach the vector coupling is approaching its 
critical value at much larger r > 0.6 fm (see Figs. I2|3i) . Our present result agrees with the 
lattice study of the function c(r) in Ref. |3l|], which will be discussed in Sect. [V] 

From Eq. (|T7|) it is evident that the asymptotic coupling a cr ; t is fully determined by the 



(21) 



ratio rf = and in two-loop approximation is given by 



a B (q = 0) 



47T 



Pi lnt c 
Po to 



with the logarithm 



t = In rj 2 = In 



Ml 



(22) 



The critical value has to be determined with great accuracy. However, small uncertainties 
in the values of Ay(n/ = 3) and Mb do not allow to fix a cr \t(nf = 3). As we have already 
mentioned, such a difference may reach ~ 30%. 
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FIG. 2: Compound «B( r ) with a cr n = 0.63015 and the parameters from Eq. (I18p (solid line), and 
ctB{ r ) (dashed line) with a cr n = 0.819 and the parameters from Eq. (|19p . 

It is important that the critical couplings in the momentum and the coordinate spaces 
coincide: 

«B(crit) = a-Q^r — > oo) = a B (? = 0). (23) 

It is of interest to understand why in phenomenological models the smaller values Ay{nj — 
3) ~ 330 — 380 MeV are often used (compare to those from (TT71)). giving, nevertheless, a 
good description of the low-lying meson states. Such values of Ay correspond to a smaller 
Ky^irif = 3) ~ 250 MeV, as compared to that from Eq. (0J, and close to those, calculated in 
the quenched approximation on the lattice, A(nf = 0) = (245 ±20) MeV 25j . Nevertheless, 
in this case a reasonable agreement with experiment is reached due to the small value taken 
for Mb and the large value of the freezing constant, a B (crit) ~ 0.60. 



IV. THE VECTOR COUPLING IN THE COORDINATE SPACE 

The vector coupling in the coordinate space is defined according to Eq. ([3]), where the 
integral can be rewritten in a different way, introducing the variable y — 

oo 

a B (rA v ,V 2 ) = - [ dy S ^^ a B (y , r, 2 ) , (24) 
7T J y 



This expression explicitly shows that «b(^) depends on the combination rAy(n/), if n/ is 
fixed, and also on the parameter rf = A ^ f ) 2 • 

In Fig. 2 two couplings qbW are shown for two sets of A(n/) from Eqs. ( 1T81) . ( TT9l) . where 
the critical values are equal to 0.630 (M B = 1.15 GeV) and and 0.819 (M B = 1.0 GeV). 
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r (GeV/c) 



FIG. 3: Comparison of the compound a.&{r) with parameters from Eq. (I18j) and «B( r ) with fixed 
nf = 3 and the same a cr ;t and Ay(n/ = 3) = 0.4974 GeV. 

In Fig. 3 the calculated "compound" a^ir) with Mb = 1.15 GeV is compared to the 
coupling qbM, in which rif = 3 is fixed (no matching), while for both couplings their 
critical values coincide and are equal to 0.630. As seen from Fig. 3, both curves are very 
close to each other and, perhaps just owing to this fact, the vector coupling with fixed rif 
may be used in phenomenological models. It also indicates that the frozen value of the 
coupling is of primary importance. 

Notice that the situation is different for light and strange mesons, which have large sizes, 
and for them a screening of the GE interaction is possible, which can occur owing to open 
channels, decreasing the vector coupling. 

V. THE STATIC FORCE AND THE FUNCTION c(r) 

To have an additional test of the calculated vector coupling «b(^) we consider here the 
static force, 

Mr) = V: t (r) =* + V' E (r) =a + (25) 

where the coupling 

a F (r) = a B {r) - ra B (r) (26) 

is introduced. The values of ap(r) are smaller as compared to a#(r) since the derivative 
a' B (r) is positive. 

In Fig. 4 the coupling «f(0 together with «b(^) with parameters from (fl~8l) are shown; 
one can see that a-p(r) is smaller by ~ 20% in the region (0.5 — 0.6) fm. 

To compare our results with the existing lattice data we introduce the dimensionless 
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FIG. 4: The compound ap(r) and ckb (r) , taking Av(rt/) from Eq. ([18]) with Mb = 1.15 GeV, 
A v (n/ = 3) = 0.4974 GeV, and a Bc rit = 0.63015. 



function r 2 FB(r) and calculate two characteristic scales: r\ and ro 32] 



where the function 



r{F n {r x ) = 1.0, r 2 F B (r ) = 1.65, 



r 2 F B (r) = r 2 cr + -a F (r), 



(27) 



(28) 



depends both on o and ap(r). For the static potential, like the Cornell and some lattice 
potentials, with the coupling constant equal to a constant, one has V st (r) = or — - (where 
e = |a(lat) = constant) and therefore in the static force, 



F(r) = o + -, 



(29) 



the coupling ctp = |e is also constant. 

On the contrary, in our calculations the coupling ap(r) changes rapidly in the region 
< r < 0.4 fm, approaching «bW only at large distances r > 0.8 fm (see Fig. HJ where 
Ay(n/) is taken from Eq. flEJ) with M B = 1.15 GeV). 

From Eq. ( 1291) one easily calculates the characteristic sizes: 



ri = 1.530 GeV -1 = 0.303 fm, r = 2.321 GeV 



-i 



0.460 fm. 



(30) 



These scales r±, r$ for Av(^i/ = 3) = 497.4 MeV, which correspond to the central value 
of the perturbative Aj^g(ri/ = 3) = 337 MeV, appear to be very close to those calculated 
on the lattice: ri(lat) = 1.621 GeV" 1 



0.321 fm 



33| and r = 1.3656(20) GeV 



-i 



(0.468 ± 0.004) fm 24|, being only several percents smaller. However, to reach precise 
agreement with the lattice scales one needs to take a bit smaller A^jg(nj), namely, those 
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values which correspond to the lower bounds in pQCD (j3J): A^g{nf = 5, 2— loop) = 208 MeV, 
A m (n f = 4) = 279.2 MeV, A m (n f = 3) = 322 MeV. For this choice we have obtained 
Qicrit = 0-5712, n = 0.312 fm, and r$ = 0.470 fm, which coincide with the lattice scales from 



Refs. 33], 34j with an accuracy better than 3%. 



In the cases considered we have found the following values for the product r$ A^{nj = 3): 

r A m (n f = 3) = 0.764, for r = 0.470 fm, a crit = 0.5712, 

r A m (n f = 3) = 0.782, for r = 0.460 fm, a crit = 0.630, (31) 



in good agreement with the lattice results from Refs. [33|, |34| . Thus we conclude that for 
large A v (% = 3) ~ 500 MeV, with the corresponding Aj^g = 0.339 MeV, one needs to take 
a relatively large IR regulator Mb = 1.15 GeV to obtain the scales r\, tq in agreement with 
the lattice results. 

For a smaller regulator, e.g. Mb = 1.0 GeV, the scales r l5 r turn out to be smaller: 
r\ = 0.292 fm and r = 0.442 fm and the product r A^g(nf = 3) = 0.757. 

On the contrary, the large value = 0.50 fm may be obtained in two ways: either with 
larger regulator Mb > 1.30 GeV, if the "perturbative" Aj^g(n/ = 3) from Eq. (jlj) are used, 
or taking the significantly smaller value of Aj^g(?7,/ = 3) ~ 245 MeV, as in the quenched 



lattice calculation 



25) 



Thus we can conclude that the scale r cannot be considered as a universal parameter 
but depends on A and M B used. Notice, that the force F B (r) depends also on the string 
tension and our calculations were performed with o = 0.18 GeV 2 , having obtained for the 
scales T\ and tq a good agreement with the lattice results. 

An additional and very important test of the vector coupling comes from the study of 
the function c(r), defined via the second derivative of the static potential and therefore it 
does not depend on the string tension. It is determined by an(r) and its first and second 
derivatives: 

c(r) = \r" V£(r) = ~a F (r) - ^og(r)y. (32) 

The second derivative a B (r) is negative and therefore the magnitude of c(r) appears to be 
smaller than that of g«p(r). Moreover, the slope of c(r) depends on the IR regulator used. 
The behavior of c(r) is shown in Fig. [5] together with the points taken from the lattice 
calculations with rif = 2 from Ref. [35]. One can see the good agreement with both results. 
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FIG. 5: The function c(r) for Mb = 1.15 GeV and Ay(n/) from Eq. {IE]); the points with the 



errors are taken from Ref. 



311 ] . where the function c(r) was calculated on the lattice with nf = 2. 



VI. THE BOTTOMONIUM SPECTRUM AS A TEST OF a crit 

In bottomonium the centroid mass M(nl) for a given multiplet nl just coincides with the 
e.v. M(nl) of the Hamiltonian (jSJ): 



2Jp 2 + mf + y st (r) 



y? nZ = M(nl)ip nl . (33) 



Here we perform calculations with the use of the SSE, since in the relativistic case the 
accuracy is better than in the nonrelativistic and so-called einbein approximations, although 
the difference between them is only several MeV for low-lying masses M(nl) and their mass 
splittings. 

The bottomonium spectrum was calculated taking the compound a^ir) with the values 
of Av(ti/) from (]18p and also in the case with fixed rif = 3 and the same Ay (ft/ = 3) = 
497.4 MeV, but without matching. It appears that the differences between the masses 
calculated in these two cases are very small, ~ 2 — 6 MeV. Therefore it is of special 
importance to consider a different Ay- (ft/ = 3), or a different Ay^(nj = 3). Our calculations 
show that there are several mass splittings, which are most sensitive to the choice of the 



ratio r? 2 = , determining the frozen value of «b(^)- Their experimental values are taken 



from 



15 



36 



231 : 



M(2P)-M(1P) = (360.0 ±1.7) MeV, 
M(3P) - M(2P) = (280 ± 14) MeV, 

M(ID) - M(IP) = A = (264 ± 2) MeV. (34) 
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Here the centroid mass M (xb{3P] 
masses measured by the ATLAS 



) ^ (10540 ± 5) MeV is estimated from two experimental 



36] and the DO [37] collaborations. 
In Ref. 38[ it was already demonstrated that the fit to the bottomonium splittings appears 
to be sensitive to choice of the critical coupling constant. However, in the vector couplings 
the parameters were often taken in a rather arbitrary way. In particular, for the lattice 



static potential with small a\ at (r) = const = 0.306 (n/ = 3) [29] the ID — IP, IP — IS 
splittings are 40-^-30 MeV smaller than their experimental values. 

First we determine the masses of the ID multiplet. The fine-structure splittings of this 
multiplet were calculated taking the strong coupling in the spin-orbit and tensor interaction 
as for the IP-states j^, namely, a s (//Fs) = 0.38 ± 0.02 at the scale ftps ~ 1 GeV. The 
following masses were obtained, 

M (1 3 P> 3 ) = (10270 ± 2) MeV, 
M{1 3 D 2 ) = (10264 ± 2) MeV, 

M(l 3£>i) = (10254 ± 3) MeV, (35) 

so that the calculated centroid mass is by ~ 1 MeV larger than M(l 3 P2) = (101 63 ±2) MeV, 
known from experiments, M(1D) = (10164 ± 2) MeV. Here we would like to notice that in 
our calculations the splittings of the ID multiplet ( l35i) are two times smaller than those in 



lattice calculations 



39). 



Fort 



da 



40 



re nP bottomonium multiplets their spin-averaged masses are known very accurately 



4l|: 



M (IP) = (9900.0 ± 0.6) MeV, M(2P) = (10260 ± 0.7) MeV, (36) 

and therefore the mass splitting M(2P) — M(1P) = 260(2) MeV is also known with great 
accuracy. 

As seen from Tabled the ID — IP splitting is in good agreement with experiment, close 
to the experimental value for Mb = 1.10 GeV, however, at the same time the 2P — IP 
splitting increases and one needs to reach the best agreement for both splittings. Choosing 
different sets of the parameters Mb and Ay (n = 3), we have observed that a good agreement 
with experimental splittings takes place in the cases with the freezing constant a cr i t = 
0.57±0.02 and the IR regulator M B = (1.15±0.02) GeV, while the calculated spin-averaged 
masses coincide with experiment within ±(5 — 10) MeV accuracy. Notice, that the lattice 
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TABLE I: The mass splittings in bottomonium in MeV (A m (n = 3) = 325 MeV) 



State 


Mn 

m b = 


— 115 GeV 
= 4.832 GeV 




— i in GeV 
= 4.840 GeV 


pvn 1 H 


ID - IP 




259 




261 


264 ±2 


2P-1P 




371 




376 


360 ±2 


3P-2P 




288 




294 


280 ± 14 



calculations give larger D-wave masses 39J, as compared to ours, while smaller masses 
Af(lD) were predicted in Ref. [421] . 

For the 2P — IP splitting a small deviation ~ 5 MeV from the experimental number 
can be obtained, if a smaller QCD constant Aj^g(?r, = 3) = (317 ± 5) MeV is used, while 
M B = (1.15 ± 0.02) GeV is relatively large. 

We do not give here the centroid masses of the 5-states, because for calculations of 
M(1S), M(2S) one needs to take into account the nonlinear behavior of confining potential 
at short distances, r < 0.20 fm and this fact gives rise to small negative corrections to the 
S'-wave masses. Such corrections are very small for the states with I ^ 0, since their w.f.s 
are equal to zero near the origin, while the 5- wave w.f.s have their maximum values there. 



VII. CONCLUSIONS 



We have studied the vector coupling in the momentum and coordinate spaces, introducing 
the IR regulator, M B = ^/2txo = (1.06 ± 0.11) GeV, as it is prescribed in BPTh. 

For the vector coupling in momentum space we have performed the matching procedure 
at the quark mass thresholds and calculated the "vector" constants A v (n/). It appears 
that these constants correspond to A^jg(n/), which coincide with the perturbative Aj^g(n/) 
within ±5 MeV; moreover, their values weakly depend on the regulator Mb, if it is taken 
from the range (1.0 — 1.20) GeV. 

The calculated scales r and T\ are shown not to be universal numbers and r decreases 
by 6% when Mb decreases from the value Mb = 1.15 GeV to Mb = 1.00 GeV. The ratio 
a = 1.505 ± 0.02 and the product r A m (n f = 3) = 0.77 ± 0.02 are obtained. 

The choice with M B = 1.15 GeV and A^fg(nj = 3) = 322 MeV gives the best (precise) 
agreement with the lattice scales r and T\ 
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The function c(r), which is proportional to the second derivative of the static potential 
and does not depend on the string tension, is calculated, and appears to be in agreement 
with the ALPHA Collab. predictions. 

Our analysis of the bottomonium spectrum shows that the splittings M (ID) — M (IP) and 
M(2P) — M(1P) are very sensitive to the ratio rj 2 = Ay ( n =3)2 used and the best agreement 
with experiment is reached taking i] = 2.46 ± 0.04, or a cr i t = 0.57 ± 0.02. 

The following splittings for the members of the bottomonium ID multiplet are predicted: 
M(l 3 D 3 ) - M(l 3 D 2 ) = 7(2) MeV, M(l 3 D 2 ) - M(l 3 A) = 9(3) MeV. 
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